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SECTION  1 
INTRODUCTION 

Design,  evaluation,  and  testing  of  trans-ionospheric  radio  frequency  (RF) 
communication  systems  require  high  fidelity  channel  models  and  detailed  knowledge  of 
fading  channel  statistics.  Such  models  can  be  used  to  constfuct  realizations  of  the 
received  signal  for  use  in  digital  simulations  or  hardware  channel  simulators.  During 
the  design  process,  knowledge  of  channel  fading  statistics  is  used  to  develop  power 
requirements,  size  interleavers,  and  assess  performance,  for  example. 

The  design  goal  for  a  truly  robust  trans-ionospheric  communications  system  is  to 
achieve  performance  that  is  acceptable  over  the  entire  range  of  possible  fading 
conditions  from  fast,  frequency  selective,  Rayleigh  fading  to  slow,  flat  fading  including 
the  regime  between  Rayleigh  fading  and  a  non-fading  channel.  While  considerable 
effort  has  been  expended  in  the  nuclear  effects  community  over  the  past  three  decades 
to  characterize  the  Rayleigh  fading  channel  and  to  develop  Rayleigh  fading  mitigation 
techniques,  somewhat  less  effort  has  been  directed  at  robust  design  and  performance  in 
non-Rayleigh  fading.  One  reason  for  this  is  that  if  a  system  is  properly  designed  to 
successfully  operate  over  the  full  range  of  Rayleigh  fading,  then  it  is  generally  assumed 
that  it  will  also  perform  well  in  non-Rayleigh  fading. 

There  are  situations,  however,  where  this  assumption  may  not  be  valid.  For 
example,  the  performance  of  many  systems  is  degraded  in  slow  Rayleigh  fading  where 
long,  deep  fades  can  cause  tracking  loops  to  loose  lock  on  the  received  signal.  A 
natural  question  to  ask  is  the  following:  What  happens  to  receiver  performance  in 
slightly  non-Rayleigh  fading  where  the  channel  coherence  time  may  be  longer  than 
under  Rayleigh  fading?  Does  performance  degrade  further  because  the  fades  may  be 
longer  or  does  it  improve  because  the  fades  are  generally  not  as  deep? 

Effects  of  non-Rayleigh  fading  are  also  important  in  determining  the 
performance  of  syst  ,ns  that  have  not  been  designed  to  operate  under  highly  disturbed 
ionospheric  conditions.  An  obvious  question  is:  W'ill  these  systems  perform  adequately 
under  weakly  disturbed  (either  man-made  or  naturally  occuning)  ionospheric 
conditions?  This  question  should  be  addressed  before  effort  is  spent  to  needlessly 
upgrade  communications  systems  that  may  already  perform  adequately  in  weak 
scintillation  or  before  fragile  communications  systems  unexpectedly  fail. 

The  purpose  of  this  report  and  a  companion  report  [Dana,  1992b]  is  to  extend 
existing  Defense  Nuclear  Agency  (DNA)  channel  models  to  include  the  non-Rayleigh 
fading  regime.  Ideally  these  models  will  then  cover  all  possible  fading  conditions,  from 
fast,  frequency  selective  fading  caused  by  a  highly  disturbed  ionosphere  to  naturally 
occurring  slow,  non-Rayleigh  fading.  Such  models  are  needed  for  the  design,  analysis, 
and  testing  of  existing  and  new  communications  systems. 

The  problem  with  a  general-purpose  channel  model  is  'hat  the  statistics  of  non- 
Rayleigh  fading  are  not  described  by  any  single  mathematical  expression,  as  is  the  case 
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for  Rayleigh  fading.  Thus  it  is  the  intent  of  the  companion  report.  Temporal  Statistics 
of  Non-Rayleigh  Fading,  to  demonstrate  that  Rician  statistics  provide  a  reasonable 
worst  case  channel  model  in  this  regime.  In  this  report,  temporal  statistics  (i.c.,  mean 
fade  duration  and  separation)  of  Rician  fading  are  derived.  These  statistics  arc  then 
used  to  demonstrate  that  realizations  of  sampled  Rician  fading  can  be  generated  with 
the  desired  statistics.  In  Dana  [1992b]  the  temporal  statistics  of  non-Rayleigh  fading  are 
analyzed,  and  it  is  shown  that  Rician  statistics  may  provide  a  reasonable  worst  case  for 
the  cumulative  distribution,  mean  fade  duration,  and  mean  fade  separation. 

1,1  RICIAN  STATISTICS, 

Realizations  of  the  channel  impulse  response  function  generated  with  Rician 
amplitude  statistics  [Rice,  1548]  have  been  used  for  many  years  to  evaluate  system 
performance  in  the  regime  between  full  Ravleigh  fading  and  ambient  non-fading 
conditiorts.  This  appioach  is  often  used  because  it  is  easy  to  generate  a  realization  uf 
Rician  fading  from  a  realization  of  Rayleigh  fading  by  simply  adding  a  constant 
component  to  the  complex  impulse  response  function,  appropriately  re-normalized  to 
maintain  constant  power. 

However  researchers  in  the  area  of  ionospheric  physics  (see,  for  example, 
Fremouw,  Livingston,  and  Miller  [1980];  Rino  and  Fremouw  [1973];  Rino,  Livingston, 
and  Whitney  [1976];  and  Whitney,  et  al.  [1972])  have  suggested  that  Nakagami-m, 
generalized  Gaussian,  or  log-normal  distributions  may  more  accurately  describe  the 
observed  amplitude  distribution  of  RF  scintillation  caused  by  the  ambient  ionosphere. 

None  of  these  distributions  or  the  Rician  distribution  adequately  describe  the 
observed  phase  fluctuations  of  non-Rayleigh  fading.  Indeed,  the  Nakagami-m  and  log¬ 
normal  distributions  only  describe  amplitude  fluctuations.  Often  two-component 
models,  one  for  amplitude  and  another  for  phase,  are  used  to  described  the  statistics  of 
obseiwed  trans-ionospheric  signals  (see,  for  example,  Wittwer  [1980]).  However,  such 
two-component  models  may  not  accurately  reproduce  observed  amplitude-phase 
correlation  of  non-Rayleigh  fading  [Fremouw,  Livingston,  and  Miller,  1980], 

Two  important  points  should  be  noted  about  proper  design  of  robust  trans- 
ionospheric  communications  systems.  First,  performance  should  be  insensitive  to  the 
random  phase  fli'^fuations  encountered  on  a  the  link.  Second,  it  is  important  that  the 
performance  is  insensitive  to  the  differences  between  the  various  fading  distributions. 
All  arc  possible,  so  the  system  should  be  designed  to  perform  against  the  reasonable 
worst  case.  If  phase  fluctuations  are  important,  then  separate  Total  Electron  Content 
(TEC)  dynamics  models  are  available  to  stress  the  system  [Wimver,  1980;  Frasier, 
1988]. 


Recently  De  Raad  and  Grover  [1990]  undertook  a  theoretical  study  of  the 
amplitude  statistics  of  non-Rayleigh  fading  for  a  wide  range  of  ionospheric  conditions. 
They  conclude  that  (1)  none  of  these  models  is  reliable  in  general;  (2)  the  actual 
amplitude  distribution  has  a  su’ong  dependence  on  the  power  spectrum  of  the  scattering 
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ionospheric  stnunure  as  well  as  the  Fresnel  length;  and  (3)  Rician  amplitude  statistics 
provide  a  useful  "worst  case”  description  of  the  occurrence  of  deep  fades. 


Multiple  phase  screen  (MPS)  techniques  <e.g.,  Knepp  [1983];  De  Raad  and 
Grover  [1990])  can  be  used  to  generate  realizations  of  the  channel  impulse  response 
function  that  represent  direct  solutions  to  Maxwell’s  equations.  These  higher  fidelity 
realizations  exhibit  a  large  range  of  amplitude  and  phase  fluctuations  under  non- 
Rayleigh  fading  conditions.  However  De  Raad  and  Grover  [1990]  correctly  observe 
that  the  uncertainty  in  the  validity  of  MPS  lealizatiors  has  been  shifted  from  the 
amplituc;  ^nd  phase  distributions  to  the  statistical  description  of  the  scattering  medium. 

Still,  the  temptation  persists  to  use  Rician  fading  realizations  for  non-Rayleigh 
fading.  They  are  easy  to  generate  from  Rayleigh  fading  realizations,  and  they  contain 
phase  fluctuations  (albeit  fluctuations  that  differ  significantly  L.,m  observations).  By 
comparison  with  MPS  realizations  of  non-Rayleigh  fading,  De  Raad  and  Grover  [1990] 
show  that  Rician  amplitude  statistics  represc  t  a  reasonable  worst  case  for  the  observed 
cumulative  distribution  of  fades,  and  Dana  [1992b]  shows  that  Rician  temporal  statistics 
also  represent  a  reasonable  worst  case  for  the  observed  mean  fade  duration  and 
separation. 


The  purpose  of  this  report  is  to  provide  further  information  on  the  temporal 
stadstics  (mean  fade  duration  and  separation)  of  Rician  fading  and  to  define  sampling 
requirements  of  Rician  realizations  of  the  channel  impulse  response  function.  These 
analytic  results  arc  corapareu  to  measured  valUvS  fiom  a  representative  set  of  MPS 
realizations  in  Dana  [1992b]  where  utility  of  Rician  temporal  statistics  in  bounding  the 
observed  range  of  fade  durations  and  separations  in  MPS  realizations  is  demonstrated. 

1.2  SAx\lPLING  STATISTICS. 


The  original  version  of  this  report  [Dana  1988]  was  intended  to  address  three 
questions  that  arise  during  simulation  or  hiidware  testing  activities  of  communications 
links  under  Rayleigh  fading  conditions:  (1)  How  many  deccr-elaiion  times  (Tq)  per 
realization  of  the  channel  impulse  response  function  axe  necessary?  (2)  How  many 
samples  per  decorrelation  time  are  necessary?  (3)  How  should  interpolation  be  done 
between  samples?  This  report  re-addresses  these  questions  for  the  more  general  case  of 
Rician  fading,  and  more  completely  addres';es  an  additional  question:  (4)  What  is  the 
expected  variation  in  measured  parameters  a  realization? 

I  The  fourth  question  can  arise  in  at  least  two  situations.  The  efficacy  of  a 

I  realization  of  the  channel  impulse  response  function  may  be  in  question,  or  it  may  be 

necessary  to  validate  a  realization  for  use  in  hardwt  .'o  testing.  An  approach  used  by  the 
author  to  validate  realizations  is  to  measure  key  pai.' meters,  such  as  mean  power, 
,  amplitude  moments,  decorrelation  time,  and  »?un;'  of  samples  per  decorrelation 

time.  These  measured  parameters  should  agree  •  ?  ensemble  values  to  within  some 

tolerance.  The  question  is:  What  tolerance?  Dana  [1991]  partially  addresses  this 
question.  This  report  incorporates  some  recent  findings  on  the  expected  tolerance  of 
;  measured  realization  para.aieters. 

i 
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The  first  three  questions  are  answered  in  part  in  the  DNA  signal  specification 
for  nuclear  scintillation  [Wiitwer,  1980]  which  requires  a  minimum  of  100 
decorrelation  times  per  realization  and  10  samples  per  decorrelation  time.  However 
considerable  statistical  variation  in  receiver  performance  is  seen  when  the  minimum 
realization  length  is  used.  This  is  particularly  true  of  links  that  have  large  power 
margins  and  are  susceptible  to  only  the  deepest  fades.  Of  course  the  best  answer  to 
these  questions  is  to  measure  link  performance  with  realizations  of  increasing  length 
and  resolution  until  the  statistical  variation  in  the  results  from  one  realization  to  the 
next  is  acceptable.  Unfortunately  the  luxur>’  of  doing  this  analysis  ordinarily  does  not 
exist. 


The  next  higher  level  of  analysis  of  these  questions  is  to  look  at  the  statistics  of 
the  realizations.  This  is  the  approach  that  will  be  taken  in  this  repon.  The  first  order 
statistics  of  realizations  are  measured  by  calculating  amplitude  moments  and  the 
cumulative  distribution  and  comparing  these  to  ensemble  values  for  Rayleigh  fading. 
The  second  order  statistics  of  the  realizations  are  measured  by  calculating  the  mean 
duration  and  separation  of  fades. 

In  general,  the  received  signal  may  be  written  as  the  convolution  of  the  channel 
impulse  response  function  h(t,t)  with  the  transmitted  modulation  m(r): 

oo 

u(t}  =  J  h{t,r)  mit-r)  dr  .  (1.1) 

0 

In  either  software  link  simulations  or  in  hardware  channel  simulato’”,  Equation  1.1 
can  be  implemented  as  a  tapped  delay  line: 

u{t)  =  2^  h{tJAx)m{i-jAf) At  (1.2) 

7=0 


where  Af  is  number  of  taps  on  the  delay  line;  At  is  the  delay  spacing  of  the  delay  line; 
and  h(tJAt)  is  the  time  varying  complex  weight  of  the 7'*  tap.  In  a  software  simulation 
of  link  performance  lime  will  also  be  discretely  sampled  (i.e.,  t  =  kAt). 

Under  Rayleigh  fading  conditions  h(r,T)  is  a  complex,  zero  mean,  normally 
distributed  random  variable  and  thus  has  a  Rayleigh  amplitude  distribution.  It  then 
follows  from  Equation  1.2  that  u(.t)  is  also  a  complex,  zero  mean,  normally  distributed 
random  variable  with  a  Rayleigh  amplitude  distribution. 

A  complete  analysis  of  the  first  three  questions  would  consider  the  sampling 
requirements  for  each  delay  of  the  discrete  impulse  response  function  hikAtJAt). 
However  this  is  beyond  the  scope  of  this  report.  Therefore  sampling  requirements  on 
the  flat  fading  impulse  response  function  li(kAt),  where 
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will  be  addressed  in  this  report.  The  sampling  requirements  for  hikAt)  will  give  some 
indication  of  the  sampling  requirements  for  the  fre^  ncy  selective  impulse  response 
function  h(kAtJAT).  Perhaps  this  shouid  be  state^  another  way;  Sampling  that  is 
inadequate  for  h(kAt)  will  surely  U  inadequate  for  h(kAtJAT).  Thus  it  is  the  intent  of 
this  report  to  define  adequate  sampling  for  h{kAt)  and  to  infer  adequate  sampling 
requirements  for  each  delay  of  hikAt JAv). 
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SECTION  2 

TEMPORAL  STATISTICS  OF  RICIAN  FADING 


This  section  is  a  generalization  of  well-known  results  from  the  classical  work  of 
Rice  [1948,  1954,  1958]  on  the  first  and  second  order  statistics  of  Rayleigh  fading.  To 
the  author’s  knowledge,  the  extension  of  Rice’s  results  on  temporal  statistics  to  non- 
Rayleigh  fading  is  nv*.w. 

2.1  FIRST  ORDER  STATISTICS. 

Under  strong  scattering  conditions,  the  electric  field  incident  on  the  plane  of  the 
receiver  is  the  summation  of  many  '’  aves  propagating  in  slightly  different  directions 
about  the  line-of-sight.  Under  the  central  limit  theorem  of  statistics,  the  two  orthogonal 
components  of  the  electric  field  must  then  be  zero-mean,  normally  distributed  random 
variables.  It  is  assumed  that  the  two  orthogonal  components  aie  also  independent,  The 
complex  narrow-band  envelope  of  the  elecuic  field  undergoing  Rayleigh  fading  may 
be  then  represented  as 


E(t)  =  x(t)  +  /  y(0 

• 

where  x  and  y  are  independent  and  noiinally  distributed  with  zero  mean  and  standard 
deviation  <t.  The  carrier  frequency  term,  exp  (icor),  has  been  neglected  in  this 
expression.  Thus  E(t)  may  be  thought  of  as  the  output  voltage  of  a  down-converter 
where  x(t)  is  the  in-phase  component  and  y(t)  is  the  quadrature-phase  component. 

Under  mild  to  weak  scattering  conditions,  a  model  of  the  electric  field  is  a 
specular  component  plus  a  normally  distributed  random  component.  The  elecuic  field 
is  then  written  ar 


£(r)  =  (x(r)  +  r  cos  t?J  +  i  [y(/)  +  t  sin  t?] 

where  r  is  the  constant  component  and  is  a  constant  phase.  Clearly  Rayleigh  fading 
corresponds  to  the  case  where  r  is  zero. 

Rice  [1948]  was  the  first  to  show'  that  the  probability  density  function  of  the 
amplitude  of  £(r). 


ait)  =  V  £(/)£*(/) 
has  the  probability  density  function: 


r  1 

ar 

=  -^expj 

lo 

y. 

where  Io(  )  is  the  modified  Bessel  function. 

For  the  mean  power  of  of  the  electric  field, 

Pq  =  (o^)  =  2(^  +  , 


(2.1) 
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to  be  constant,  the  power  of  the  fluctuating  component,  2cr^,  must  be  reduced  as  the 
power  of  the  specular  component,  r‘‘,  is  increased.  To  keep  track  of  this  in  a  consistent 
manner  both  are  written  in  terms  of  the  scintillation  index  54,  where 


= 


,2v2 


The  powers  of  the  two  components  are  then: 

=  Pq{\-R) 
^  PoR 

where  the  “Rician”  index  R  is 

R  =  Vl-54  . 


The  Rician  index  is  the  fraction  of  the  total  power  that  is  in  the  constant  component. 
Upon  writing  crand  r  in  terms  of  R,  Equation  2.1  becomes 

2^^Ra^/Po\ 


2a 


a^/Pg  +  R 
■  \-R 


The  corresponding  phase  6  ;f ‘he  electric  field  is 

\x(t)  +  rcos  i? 


lo 


1-/? 


0  =  tan 


-i 


y(0  +  rsin  1? 


The  probability  density  function  of  the  phase  is 


1 

R 

1 

R 

-  2k 

"l-R 

L.  J 

^2 

k{\-R)_ 

tL 

2 


cos(0-t?)  X 


P[1  -cos(0-i?)  ] 
^^p]-  - iZr 


1  +  erf 


'  R 

LV 


l-R 


cos(0- 1>) 


where  erf  (•)  is  the  error  function.  The  probabilTy  density  function  of  phase  is  just 
l/2;i  when  the  scintillation  index  is  unity  (/?=0),  as  it  should  be  for  Rayleigh  fading. 

These  two  probability  density  functions  are  plotted  in  Figures  1  and  2  for 
several  values  of  the  scintillation  index.  Amplitude  in  Figure  1  is  oAfp^,  and  phase  in 
Figure  2  is  the  quantity  0-i?.  The  solid  line  curves  in  both  figures  are  the  Rayleigh 
limits.  As  expected,  both  functions  approach  delta  functions  as  the  scintillation  index 
approaches  zero. 
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Figure  1.  Rician  amplitude  probability  density  fmiction 
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f  igure  2.  Rician  phase  probabiiity  density  function. 
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The  cumulative  distribution  of  the  power  P  {P  =a^),  which  is  equal  to  the 
probability  that  the  instantaneous  power  is  less  than  or  equal  to  P,  is  given  by 

r  1 

F(P)  =  ]f{a)da  =  exp  77^  j  j ^\jZr  J 

where  y{n,x)  is  the  incomplete  gamma  function  and  r(/i+l)  =  n!  is  the  gamma 
function.  The  summation  is  obtained  by  expanding  the  Bessel  function  in  a  power 
series  and  then  performing  the  integration  term-by-term.  This  form  of  Marcum’s  Q 
function  [Marcum,  1948]  is  easily  evaluated  for  values  of  R  that  are  not  too  close  to 
unity.  In  particular  Equation  2.2  converges  slowly  for  values  of  54  less  than  0.25, 
corresponding  to  R  values  greater  than  0.96. 

For  the  Rayleigh  case  the  cumulative  distribution  is  exponential: 

F(/>)  =  1-exp  .  (54=1) 

.  ^0  J 

The  Rician  cumulative  distribution  function  is  plotted  in  Figure  3  versus  the 
ratio  P/Pq  for  several  values  of  the  scintillation  index.  For  values  of  the  scintillation 
index  between  0.75  and  unity  the  Rician  cumulative  distribution  is  close  to  the  Rayleigh 
curve.  As  the  scintillation  index  is  reduced  from  about  0.75  to  0.5,  the  probability  of 
deep  fades  is  sigoiftcantiy  reduced.  It  is  noteworthy  that  case  where  the  power  of  the 
specular  and  fluemating  components  are  equal  corresponds  to  an  54  value  of  V374  (54  = 
0.866).  Thus  a  Rician  cumulative  distribution  does  not  deviate  significantly  from  a 
Rayleigh  distribution  until  more  than  half  of  the  power  is  in  the  constant  component. 
Of  course  the  performance  a  receiver  may  be  quite  sensitive  to  the  existence  ol  a 
specular  component. 

2.2  DBPSK  EXAMPLE. 

An  easily  calculated  example  of  the  effects  of  Rician  fading  is  the  differentially 
coherent  binary  phase-shift  keying  (DBPSK)  symbol  error  rate.  The  well  known 
DBPSK  symbol  error  rate  for  an  additive  white  Gaussian  noise  (AWGN)  channel  is 

PsE  =  \e-yP  (AWGN  Channel)  (2.3) 

where  yis  the  symbol  energy-to-noise  density  ratio  and  P  is  unity  for  this  cnannel.  In  a 
Rician  fading  channel  this  error  rate  must  be  averaged  over  the  probability  density 
function  of  the  fading  power  P  -  a^: 

00 

r  I  7 

(Pse)  =  ■ 

J 

0 
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Figure  3.  Cumulative  distribution  of  Rician  fading. 


This  equation  is  easily  evaluated  using  Equation  2.1  with  the  result: 

1  |”  ~\ 

{Pse)  =  2(7+ (T^H  (Rician  Channel)  . 

When  R  is  unity,  corresponding  to  the  non-fading  case,  this  expression  reduces  to 
Equation  2.3.  When  R  is  zero,  corresponding  to  full  Rayleigh  fading,  it  induces  to  the 
well-known  form: 


{Pse)  =  2[i +Y]  (Rayleigh  Channel)  . 

Plots  of  the  Rician  channel  DBP.SK  error  rates  for  several  values  of  the 
scintillation  index  are  in  Figure  4.  As  one  might  expect  from  examining  the  cumulative 
distribution,  the  DBFSK  symbol  eaor  rate  for  a  Rician  fading  channel  is  close  to  the 
full  Rayleigh  fading  channel  error  rate  when  the  scintillation  index  is  larger  that  about 
0.75,  and  is  close  to  the  AW'GN  error  r.ate  when  the  scintillation  index  is  less  that  about 
0.25.  Thus  for  DBPSK  the  most  interesting  values  of  scintillation  index,  excluding  1 
and  0,  are  between  0.75  and  0.25. 
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Symbol  Energy-to-Noise  Density  Ratio,  y(dB) 

Figure  4.  DBPSK  symbol  error  rate  for  Rician  fading. 


2.3  SECOND  ORDER  STATISTICS. 

The  fading  rate  is  determined  by  the  second  order  statistics  of  the  fluctuating 
p  it  of  the  electric  field.  The  autocovarience  of  the  electric  field  is,  in  general, 

<[£:(/)-£o][^Vt)-4]>  =  <.r(r)4/+T)) +<y(r)y(/+r)>  =  2<T^p(T) 


where 


f  0  -  r  cos  t?  +  /  r  sin  t?  . 


There  are  two  limiting  forms  for  the  correlation  function  p(T).  Under  strongly 
disturbed  scattering  conditions  that  occur  at  early  times  or  ai  the  center  of  the 
disturbed  region,  p(r)  has  the  Gaussian  form 


p(T)  =  exp 


11 


where  Tq,  the  decorrelaiion  time  of  the  electric  field,  is  defined  as  the  e  folding  point 
of  the  autocorrelation  function  [p(To)  =  e"*]-  The  corresponding  Doppler  spectrum  of 
the  temporal  fluctuations  is 


S{ojd)  = 


cc  _ 

J  exp  (-itOpT)  p(T)  dz  -  ^fnzo  exp 

— oo 


4 


which  also  has  the  Gaussian  form.  Under  less  disturbed  conditions,  the  correlation 
function  is  usually  assumed  to  have  the  form 


p(T)  = 


0411! 
^0  . 


where  the  parameter  Ui  =  2.146193)  is  determined  by  the  condition  that  p(To)=e“^ 
The  concpponding  Doppler  spectrum  has  the  form  commonly  referred  to  as  an  /  ^ 
spectrum: 


04  [l  +  itocop/a^ff 

A  third  Doppler  spectrum  is  used  for  real-time  frequency  selective  channel 
models  [Dana  1992a].  This  /'*  spectrum  has  the  functional  form 

^  16Tft _ 1 _ 

y.coo)  =  (tocoofaef? 

where  *he  normalization  of  S(<Uo)  is  chosen  so  that  p(0)  is  unity.  The  corresponding 
correlation  function  is 


P(T)  = 


a6lTl  (a^zf 

^0  3to 

exp 

'To  . 

where  =  2.904630  results  from  setting  p(To)=e  ^ 


A  comparison  of  Rayleigh  fading  realizations  (54=  1)  of  the  impulse  response 
function  with  and  Gaussian  Doppler  spectra  is  shown  in  Figure  5  where 

realization  power  in  decibels  (dB)  is  plotted  versus  time/To.  These  realizations  were 
generated  from  the  same  set  of  random  numbers,  as  described  in  Appendix  A,  so  there 
is  correlation  in  the  features  seen  in  the  three  frames.  The  realization  in  the  bottom 
frame  has  the  most  spiky  appearance  because  it  has  more  energy  at  high  Doppler 
frequencies.  The  three  realizations  have  similar  low  frequency  behavior,  and  fades  in 
the  realizations  follow  each  other  quite  closely.  The  difference  between  the  realizations 
is  the  high  frequency  jitter  of  the  and  f~^  realizations  about  the  more  smoothly 
varying  Gaussian  one.  The  significance  of  this  on  the  temporal  statistics  of  the  fades 
will  become  apparent  later. 

I 
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A  comparison  of  Rician  fading  realizations  of  the  impulse  response  function  with 
Doppler  spectra  and  scintillation  indices  of  1.0,  0.75,  0.5  and  0.25  is  shown  in 
Figure  6.  Again  these  realizations  were  generated  from  the  same  set  of  random 
numbers,  so  there  is  correlation  in  the  features  seen  in  all  four  plots.  As  expected,  the 
Rayleigh  fading  realization  in  the  top  frame  (reproduced  from  the  bottom  frame  of 
Figure  3)  has  the  deepest  fades  and  the  largest  flares.  As  the  scintillation  index  is 
reduced,  the  deep  fades  fill  in  and  the  power  in  flares  above  0  dB  is  reduced.  It  is 
interesting  that  as  the  scinti^ation  index  decreases,  fades  at  a  given  level  appear  to  get 
longer.  This  phenomenon  is  shown  theoretically  for  Rician  fading  in  the  developments 
below. 

2.4  TEMPORAL  STATISTICS. 

The  mean  duration  and  separation  of  fades  below  an  aibitrary  power  level  P  and 
that  of  flares  above  P,  are  calculated  from  the  mean  number  {NiP,T))  of  crossings  of 
the  level  P  in  the  time  interval  T. 

The  probability  that  the  amplitude  a  crosses  the  level  /  =  Vp  in  the  time  interval  t 
to  t-¥dt  with  a  positive  derivative  is  equal  to  the  probability  that  a'  >  0  and  that  I  -  a'dt 
<a<l.  This  probability  is  given  by  the  expression 

OO  I  oo 

J  da'  f  dajia,a')  =  dt  j  da' a' J{1, a) 

0  /-fl'dr  0 

where  f{a,a')  is  the  joint  probability  density  function  of  the  amplitude  a  and  its  time 
derivative  a'  =  da/dt.  The  probability  that  a  will  cross  the  level  /  in  the  time  interval  t 
to  (+dt  with  a  derivative  of  either  sign  is  then 

oo 

dt  J  \al  f(/.«')  da'  . 


For  stationary  processes,  the  mean  number  of  level  crossings  of  P  in  the  interval  t  to 
t+T  then  becomes 


OO 

{N{P,T))  =  T  J  la'I  f(VP,fl')  da'  . 

~-oo 


The  joint  probability  density  function  of  the  Rician  distributed  amplitude  a  and 
its  time  derivative  a'  is  derived  in  Appendix  B.  This  function  is: 
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Figure  6.  Realizations  of  Rician  fading  with  Doppler  spectra  and 
scintillation  indices  of  1.0,  0.75,  0.5  and  0.25. 
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It  can  be  seen  from  the  form  of  this  equation  that  the  probability  density  function  of  a 
is  Rician;  the  probability  density  function  of  a'  is  Gaussian  with  zero  mean  and 
variance  of  and  a  and  a'  are  independent.  Also  the  functional  form  of 

f{a,a)  is  independent  of  the  functional  form  of  the  Doppler  spectrum.  Only  the 
parameter  4  varies  with  the  Doppler  spectrum  (4  =  1  for  the  Gaussian  spectrum,  4  = 
1.1858  for  the/*^  spectrum,  and  4  =  1.518  for  the  spectrum). 

The  mean  number  of  level  crossings  can  now  be  easily  evaluated  with  the  result 


(Af(P,T)>  =  4 


r  8p/pq  1 

L 

2 

'  P/PQ+P' 

Jo, 

,It(i-p)J 

oxp 

■  1-P 

lo 


2yjRP/Pa 

l-R 


The  effect  of  different  Dopnler  spectra  is  to  scale  the  mean  number  of  level  crossings 
by  the  quantity  4.  This  fact  was  shown  qualitatively  by  comparing  the  realizations  with 
different  snectid  in  Figure  5. 

Figure  7  shows  plots  of  the  mean  number  of  crossings  of  P  in  one  decorrelation 
time  versus  the  ratio  P/Pq  for  a  Gaussian  Doppler  spectrum  and  several  values  of  the 
scintillation  index.  For  the  Rayleigh  case  the  maximum  value  of  <//(P,D>  occurs  at 
P/Pq  =  1/2  or  -?  dB.  As  the  scintillation  index  decreases  the  maximum  value  of 
<N(P,T))  approaches  0  dB. 


By  noting  that  two  level  crossings  are  required  to  define  the  beginning  and  end 
of  a  fade,  the  number  of  fades  per  unit  time  below  the  level  P  is  T]  =  {N{P, Tq))/2Tq. 
The  mean  separation  {Tsep{P))  of  fades  below  P  is  then  obtained  from  the  mean 
number  of  fades  per  unit  time.  For  any  long  time  interval  T  the  mean  number  of  fades 
is  TfT,  and  the  mean  separation  is  just  T/riT  or  l/rj.  Thus  the  mean  separation  of  fades 
below  P  is 


(Ts.piP))  -  • 

The  mean  separation  of  fades  below  P  is  equal  to  the  average  time  between  crossings  of 
P  with  either  a  negative  value  of  a'  (which  defines  the  start  of  the  fade)  or  with  positive 
value  of  a'  (which  defines  the  end  of  the  fade).  Thus  the  mean  separation  of  fades 
below  P  is  also  equal  to  the  mean  separation  of  flares  above  P. 
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Figure  7.  Mean  number  of  level  crossings  per  Tq  for  several  scintillation 
indices. 


The  mean  duration  (JouAP))  of  fades  below  P  is  obtained  as  follows:  During  a 
long  time  interval  T  the  total  time  that  the  power  will  be  below  P  is  F{P)T  where  F(,P) 
is  the  cumulative  distribution  given  in  Equation  2.2.  The  mean  duration  is  then  the  sum 
of  all  durations  F{P)T  divided  by  the  number  of  fades  t]T.  The  result  is 

{lDur{P))  -  ^SiP,To))  ■ 


The  mean  duration  {TfianiP))  of  a  flare  above  P  is  the  mean  time  that  the  power 
stays  above  P.  Using  the  arguments  given  above,  the  mean  separation  of  a  fade  or  a 
f.are  is  equal  to  the  mean  time  that  the  signal  is  above  P  plus  the  mean  ume  that  it  is 
below  P,  (TouriP))  +  (TfiareiP))  =  (TsepiP))-  The  mean  duration  of  a  flare  is  then 


{TflarciP))  = 


2tb[l-F(P)] 

<NiP,ro)} 


The  mean  duration  and  separation  of  fades  are  shown  in  Figures  8  and  9, 
respectively,  for  a  Gaussian  Doppler  spectrum  and  several  values  of  the  scintillation 
index.  For  other  Doppler  spectra,  the  curx'es  in  Figures  8  and  9  scale  by  IM. 

The  curves  in  Figure  8  show,  for  some  power  levels,  that  the  durat'on  of  fades 
increases  as  the  scintillation  index  is  reduced.  The  mean  duration  of  fades  for  equal 
to  0.75  exceeds  that  of  Rayleigh  fading  at  all  power  levels;  the  mean  fade  duration  for 
Si  equal  to  0.5  exceeds  that  of  Rayleigh  lading  except  for  P/Pq  values  between  -4  and 
-2  dB;  and  the  mean  fade  duration  for  54  equal  to  0.25  exceeds  that  of  Rayleigh  fading 
except  for  P/Pq  values  between  -13  and  -1  dB.  Note,  however,  that  when  is  equal  to 
0.25  the  probability  of  a  13  dB  fade  is  3.7x10'*^  and  the  mean  separation  of  13  dB 
fades  is  7.4x1  O^Tq.  Thus  for  ^4  values  greater  than  about  0.25,  it  is  possible  to  have 
fades  that  are  longe;  than  occur  at  the  same  level  with  full  Rayleigh  amplitude 
statistics. 


*35  -30  -25  -20  -15  -10  -5  0  5  10 

Power  Level,  P/Pg  (dB) 


Figure  8.  Mean  duration  of  Rician  distributed  fades. 
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SECTION  3 

SAMPLED  RiCIAN  FADING 


The  requirements  on  the  sampling  of  fading  realizations  are  given  in  the  DNA 
signal  specification  for  nuclear  scintillation  [Wittwer,  1980]  which  requires  a  minimum 
of  100  decorrelation  times  per  realization  and  10  samples  per  dcconelation  time.  The 
questions  that  arise  from  this  requirement  can  be  summarized  as:  How  close  are  such 
realizations  to  Rician  fading?  To  address  this  question,  random  realizations  of  Rician 
fading  are  generated;  moments  of  the  amplitude,  cumulative  distribution,  and  mean 
fade  duration  and  separation  are  measured;  and  these  measured  values  are  compared 
with  their  ensemble  values. 

Dana  [1988]  showed  that  lOOib  realizations  of  Rayleigh  fading  are  adequate  for 
fade  depths  of  20  dB  or  less,  and  that  40010  realizations  are  necessary  to  simulate  fades 
down  to  30  dB.  It  was  also  shown  in  Dana  [1988]  that  10  samples  per  decorrelation 
time  are  sufficient  when  linear  interpolation  of  the  complex  impulse  response  function 
is  used  to  sample  h{kAtsam)  40  times  per  decorrelation  time. 

Thus  results  in  this  section  are,  for  the  most  part,  limited  to  400%  realizations 
sampled  at  Only  the Doppler  frequency  PSD  is  considered  for  non- 

Rayleigh  fading  because  this  is  the  PSD  recommended  by  DNA  for  slow,  flat  fading 
cases  where  the  scintillation  index  is  most  likely  to  be  less  than  unity. 

Because  of  the  finite  number  oi  samples  iu  each  realization,  each  measurement 
of  realization  statistics  is  a  random  variable  with  some  mean  and  standard  deviation. 
Variations  in  statistics  from  realizaiion-to-realization  are  measured  by  generating  a 
large  number  of  realizations  (1024  to  be  exact).  Each  parameter  is  measured  by 
averaging  over  the  entire  realizatioi .  Average  and  standard  deviation  values  of  the 
1024  measurements  are  computed.  Thus  the  standard  deviations  below  represent  the 
realization-to-realization  variation  in  the  measurements  of  amplitude  moments, 
cumulative  distribution,  and  temporal  statistics. 

The  measurement  variation  of  the  mean  power  of  realizations  can  be  calculated 
analytically,  as  discussed  in  Appendix  C.  It  may  be  possible  to  compute  measurement 
variances  for  other  amplitude  moments  in  the  general  case  of  Rician  fading.  Such  a 
tedious  exercise,  however,  is  left  to  the  determined  reader.  Power  measurement 
variances  below  agree  quite  well  with  the  analytic  results  given  in  Appendix  C. 

Three  cases  will  be  considered.  The  number  of  samples  per  realization  N  is 
1024,  2048,  or  4096,  and  the  number  of  samples  per  decorrelation  time  Nq  is  10.  Here 
Nq  is  the  number  of  samples  per  decorrelation  time  used  to  generate  the  realizations. 
Methods  of  generating  such  realizations  are  outlined  in  Appendix  A.  To  measure  the 
statistics  of  the  realizations,  linear  interpolation  of  the  real  and  imaginary  paits  of  the 
impulse  response  function  is  used  to  obtain  a  sampling  period  Atscu-AOf  %/40. 

The  objectives  of  this  section  are  to  present  the  means  and  standard  deviations  of 
amplitude  moment,  cumulative  distribution,  and  temporal  statistics,  and  to  attempt  to 
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answer  the  above  question  based  on  these  results.  This  section  is  limited  to  flat  fading 
realizaiions  with  varying  values  of  the  scintillation  index  and  Doppler  frequency  power 
spectral  density. 

3.1  MEASURED  FIRST  ORDER  STATISTICS. 

One  criterion  for  deciding  that  a  realization  has  the  proper  Rician  amplitude 
statistics  is  that  measured  moments  of  the  amplitude  should  agree  with  Rician  values 
within  some  tolerance.  Ensemble  values  for  the  moments  of  the  amplitude  are  obtained 
from  Equation  2.1; 


H  (2-R^) 

where  Io(  )  and  Ii(  )  are  modified  Bessel  function. 


Tnese  moments  are  plotted  in  Figure  10  versus  the  scintillation  index  for  unity 
mean  power.  However,  amplitude  moments  are  easily  obtained  for  other  values  of  tiie 
mean  power  by  noting  that  <a”>  scales  as  P'^. 

The  scintillation  index  S4  is  the  standard  deviation  of  the  power.  It  is  necessary 
but  not  sufficient  that  ^4  equal  unity  for  Rayleigh  fading.  The  scintillation  index  is  a 
good  measure  of  the  statistics  of  flares  but  not  of  fades. 

Statistics  that  are  sensitive  to  the  distribution  of  fades  are  moments  of  the  log 
amplitude.  Using  Equation  2.1  and  a  little  algebra,  these  moments  are  found  to  be: 

<x)  =<1M)  ' 
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Scintillation  Index, 

Figure  10.  Amplitude  moments  of  the  Rician  distribution. 


=(ln2<i>  =  ln[P(,(l-R)]<x)-jln^[Po(l-«)] 


The  V  and  ^  functions  a.e: 

^  1  1 

\|/(n+l)  =  -Y  +  ^  =  V(n)  +  - 

V2.n+\)  =  =  C(2.-i)-4.  ?(2,1)  = 

k=n+ 1  AC  n  o 

where  y  is  Euler’s  constant  (y  =  0.5772157  -). 
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The  first  two  moments  of  log  amplitude  are  plotted  in  Figure  1 1  versus  the 
scintillation  index  for  unity  mean  power. 

Measured  values  of  the  mean  and  standard  deviation  of  the  amplitude  moments, 
Sn,  <X>,  and  (x^>  for  Rician  fading  realizations  are  in  Table  1  for  eight  cases  including 
lOOto,  200ro.  and  400To  long  realizations,  three  different  Doppler  frequency  power 
spectra]  densities  (PSDs),  and  four  values  of  the  scintillation  index.  Measured  values 
for  a  single  realization  should  equal  the  ensemble  value  plus  or  minus  one  or  two 
standard  deviations.  It  can  be  seen  from  the  table  that  the  average  values  are  close  to 
the  ensemble  values  but  the  standard  deviations  of  the  higher  amplitude  moments  can 
be  as  large  as  20  percent  of  the  measured  values. 

It  is  noteworthy  that  the  measurement  variation  of  (%)  increases  dramatically  as 
the  scintillation  index  is  reduced  while  that  of  (x^)  is  relatively  insensitive  to  54.  An 
explanation  for  this  curious  behavior  has  not  been  discovered. 


Scintillation  Index,  S^ 


Figure  11.  Log  amplitude  moments  of  the  Rician  distribution. 


Table  1. 

Statistics  of  sampled  Rician  fading  realizations. 


Ensemble  Values 


1024 

2048 

4096 

4096 

4096 

4096 

10 

10 

10 

10 

10 

10 

Gauss 

Gauss 

Gauss 

/•^ 

1.0 

1.0 

1.0 

1.0 

1.0 

0.75 

Measured  Values*  (Normalized  to  Ensemble  Values 


0.999 

0.054 


.998 

.105 


0.991 

0.226 


r\ 

0.084 


1.000 

0.213 


.996 

.153 
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83 


0.997 

0.028 


0.994 

0.054 


0.993 

0.0117 


1.997 

0.991 

i.085 

0.081 

0.061 


1.000 

0.111 


1.011 


1.004 

0.114 


1.001 

0.081 


1.004 


0.056  0.038 


*  p  =  Measured  average  value 
o  =  Measured  standard  deviation 


0.044 


1.013 

0.109 


1.007 

0.075 


0.998 

0.048 


Perhaps  a  better  criterion  for  the  validity  first  order  statistics  is  close  agreement 
between  the  Rician  and  the  mcasiued  cumulative  distributions.  Measured  cumulative 
distributions  (dots  plus  or  minus  one-sigma  error  bars)  are  plotted  in  Figures  12-16 
for  cases  5-8  in  Table  1 ,  respectively,  along  with  the  ensemble  curves  (Eqn.  2.2).  A 
level  of  0  dB  corresponds  to  the  mean  power  Pq.  It  can  be  seen  from  the  figures  that 
400ih  realizations  do  indeed  have,  on  the  average,  a  Rician  distribution  of  fades. 
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3.2  MEASURED  SECOND  ORDER  STATISTICS. 

Table  1  also  contains  the  mean  and  standard  deviation  of  the  measured  number 
of  samples  per  decorrelation  time.  The  measured  value  of  Nq  is  obtained  by 
performing  an  autocorrelation  of  the  complex  impulse  response  function  and  finding 
the  e“‘  point.  Close  agreement  of  this  parameter  with  its  ensemble  value  ensures  that 
the  realization  will  indeed  have  the  desired  decorrelation  time  in  a  simulation  or 
hardware  test. 

The  fidelity  of  the  realizations  in  reproducing  the  second  order  statistics  of 
Rician  fading  will  be  demonstrated  by  considering  the  mean  fade  duration  and 
separation.  The  mean  fade  duration  is  a  good  statistic  to  examine  for  communications 
applications  because  errors  often  occur  in  bursts  during  deep  fades.  If  the  fades,  on  the 
average,  are  too  long  or  too  short,  error  bursts  will  not  have  the  proper  durations  and 
the  resulting  receiver  performance  may  be  misleading. 

Fade  duration  and  separation  measurements  (dots  plus  or  minus  one-sigma  error 
bars)  and  ensemble  curves  (solid  lines)  for  the  f~*  Doppler  PSD  are  shown  in  Figures 
16-19  for  cases  5-8  respectively.  Figure  16  shows  the  these  measurements  for  full 
Rayleigh  fading.  Good  agreement  between  the  measure  and  ensemble  values  is  seen, 
except  for  30  dB  fades.  At  this  level  the  ensemble  fade  duration  is  0.026  Tq,  which  is 
quite  close  to  the  sample  duration  of  0.025 Tq.  As  the  scintillation  index  is  reduced, 
measured  mean  fade  durations  are  generally  close  to.  if  not  right  on,  the  ensemble 
curves.  However,  large  variations  are  seen  in  the  measured  mean  fade  separations. 

Two  effects  contribute  to  the  low  mean  fade  separation  measurements  seen  for 
scintillation  indices  less  than  unity.  Because  separation  measurements  require  two 
fades,  some  realizations  do  not  contribute  to  the  realization-to-realization  fade 
separation  statistics,  thereby  reducing  the  measured  average.  Also,  it  is  nor  possible  in 
these  realizations  io  measure  fade  separations  larger  than  about  400To.  Thus  measured 
mean  separations  are  biased  to  lower  values  because  large  random  samples  are  absent. 

In  Figure  19,  for  a  scintillation  index  of  0.25,  a  large  discrepancy  in  the 
measured  mean  fade  duration  and  separation  is  seen  at  fade  level  of  +3  dB.  Here  the 
mean  fade  duration  is  about  SOOrh.  which  clearly  cannot  be  accurately  measured  in 
400 To  realizations.  Because  the  measured  mean  fade  separation  and  measurer  mt 
variance  at  this  level  arc  not  zero,  at  least  two  realizations  must  have  had  two  +3  dB 
flaies  even  though  the  probability  of  such  an  event  is  less  than  5x10“'*.  The  measured 
mean  duration  of  5  dB  fades  (which  occur  with  a  probability  of  approximately  3x10“'*) 
is  close  to  the  ensemble  value  of  0.3  IIq,  but  the  measured  mean  fade  separation  is  only 
ten  percent  of  the  ensemble  value  of  1042To. 

Except  for  low  probability  events  with  large  durations  or  separations,  these 
results  demonstrate  that  400To  Rician  fading  realizations  sampled  at  To/40  do  indeed 
have  the  proper  temporal  statistics. 


25 


I 

I 

I 


Figure  14.  Measured  cumulative  distribution  for  54  =  0.5. 
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Figure  15.  Measured  cumulative  distribution  for  54  >  0.25. 
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Figure  16a.  Measured  mean  fade  duration  for  54  s  1.0. 
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I  Figure  16b.  Measured  mean  fade  separation  for  54  =  1.0. 
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Figure  18a.  Measured  mean  fade  duration  for  S4  : 
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Figure  18b.  Measured  mean  fade  separation  for 
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Figure  19a.  Measured  mean  fade  duration  for  S4  =  0.25. 
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Figure  19b.  Measured  mean  fade  separation  for  ^4  =  0.25. 
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APPENDIX  A 

REALIZATIONS  OF  RICIAN  FADING 

The  methods  of  generating  flat  fading  realizations  of  the  impulse  response 
function  for  Rician  fading  channels  are  simple  extensions  of  Rayleigh  fading  methods 
discussed  in  detail  elsewhere  [e.g.,  Knepp,  1982;  Knepp  and  Wittwer,  1984\  Dana, 
1991;  Dana  1992a].  Only  a  brief  review  is  presented  this  appendix. 

The  methods  fall  into  two  classes:  Fourier  transform  and  real-time  digital  filter 
techniques.  The  f~^  and  Doppler  Power  Spectral  Density  (PSD)  realizations  are 
particularly  I'mple  to  £  erate  using  digital  filters.  Realizations  with  a  Gaussian 
Doppler  PSD  are  more  easily  generated  using  Fourier  transforms.  These  methods  will 
be  reviev/ed  in  subsections  of  this  appendix. 

A,1  FOURIER  TRANSFORM  TECHNIQUE  FOR  GAUSSIAN  PSDs. 

The  starting  point  of  the  Fourier  transform  method  of  generating  a  realization 
of  flat  Rician  fading  is  the  Doppler  frequency  PSD  function,  S{C0d).  The  Gaussian 
form  of  this  function  is: 


S{(Od)  =  exp 


4  . 


(Gaussian  PSD) 


where  Pq  is  the  mean  power  of  the  realization,  and  \-R  is  the  fraction  of  power  in  the 
random  component.  The  "Rician”  index  R  is 

R  =  ^l-Sl 


where  S4  is  the  scintillation  index.  The  quantity  SiQ)D')dcoD/2Ti  is  the  mean  power  in  the 
Doppler  radian  frequency  interval  coed^K  to  {(Oo  +  dwoV^n. 

Discrete  realizations  of  the  channel  impulse  response  function  will  contain  N 
time  samples  and  Nq  samples  per  decorrelation  time.  Thus  the  time  spacing  of  the 
discrete  samples  is 


and  the  total  time  duration  of  the  realization  is  N^t,  In  the  Doppler  radian  frequency 
domain  the  spacing  of  the  discrete  samples  is 

^OlD  =  . 

Note  that  the  quantity  Acood^t,  which  v/ill  appear  later  in  a  Fourier  transform,  is  just 
2k/N. 


A-I 


The  samples  in  the  frequency  domain  are  generated  by  first  calculating  the 
fraction  o  signal  power  in  each  Doppler  frequency  bin,  Sj  =  S{jAa)D)A(ODl2n.  For  the 
Gaussian  PSD. 


^  _  4^P^{\-R)Na  r  /nX" 

Sj  -  ^  expj^  ^2 


(/  =  -W2,  -  ,W2-1)  . 


Next  the  random  Doppler  frequency  spectrum  HijAcoo)  of  the  impulse  response 
function  is  generated: 


where  R  is  the  fraction  of  the  total  power  in  the  constant  component  and  ^  is  the 
constant  phase  of  the  Rician  component.  The  quantity  Sj,k  is  the  Kronecker  delta 
symbol: 


k 


1  J  =  k 
G  otherwise 


The  leading  factor  IkIAqjd  has  been  included  in  HQAcod)  so  that  the  discrete  Fourier 
transform  of  HijAcoo)  will  be  dimensionless.  Random  components  of  the  spectrum,  . 
are  complex,  normally  disturbed  random  variables  with  the  properties: 

(^J  ^k)  =  ^J.k 
^k)  =  0  . 

Thus  the  mean  power  of  the  ^  samples  is  unity.  Random  samples  of  may  be  easily 
generated  using 


=  V-ln(wi^-)  exp  i2Kiu2j) 

where  u^j  and  U2j  are  independent  random  variables  uniformly  distributed  on  the 
interval  [0,1). 

Finally  the  random  Doppler  spectrum  of  the  channel  impulse  response  function 
is  Fourier  transformed  to  the  time  domain.  In  continuous  notation  this  Fourier 
transform  is 


Ht) 


oo 

H{0)d)  exp  iiojot) 

» 

-*oo 


dcOp 

2jz 


and  in  discrete  notation. 


A-2 


A72-1 

V  ^COn 

h{kAt)  =  ^^^HQAcod)  exp  {i{jA(OD)i.kAt)] 

NQ^\ 

=  (V^Po/?ei<p)5y.o  ]  exp  [iKHjk/^T)] 


where  k  =  0,1,  •••  JV-1. 

A.2  DIGITAL  FILTER  TECHNIQUE  FOR AND PSDs. 

An  or  realization  can  be  generated  by  passing  white  Gaussian  noise 
through  cascaded  single-pole  filters,  as  desaibed  in  Dana  [1992a]. 

An/“^  filter  can  be  created  by  cascading  two  single-pole  filters: 

yk  ayk-i  +  hxic-i  (A.l) 

xk  =  axk-i  +  bvk-i  . 


The  coefficients  a  and  arei; 


a  ~  ejpi-a^At/^)  =  exp(-a4Wo) 

b  =  nT? 


where  =  2.146193.  Discrete  samples  of  the  additive  white  Gaussian  noise  random 
piocess  Vic  are  generated  using  the  equation: 


v* 


(l~a^)Po{\-R) 

1+0^ 


L 

2 


I*  . 


The  vjt  samples  must  have  mean  power  given  by  the  quantity  in  the  square  brackets  so 
that  the  mean  power  of  the  filter  output  samples,  yk,  will  have  mean  power  of 
The  discrete  channel  impulse  response  function  in  this  case  is 

h{kAt)  =  yk  +  ^jPoRe^V  . 

To  minimize  the  transient  response  at  start-up  it  is  necessary  to  initialize  the 
filter.  This  is  done  by  setting 


*  Some  authors  prefer  to  include  the  gain  of  tlic  filter  in  b  coefficients.  For  example,  see  Bogusch  [  i989] 
Equation  2-40.  Wittwer  [1980]  writes  the /~*  filter  equations  as  shown  here  but  he  combines  the 
exponential  in  the  a  coefficient  with  the  expression  for  Uie  mean  power  of  the  input  white  Gaussian  noise 
to  obtain  a  hyperbolic  tangent  function. 


A-3 


yo  -  '^Pq{^-R)i*q 


-<^0 


{l-o^)Po(l~R) 

IV 


where  uq  and  Uj  are  independent  samples  of  the  random  process  ^  uncorrelated  with 
the  vie  samples.  The  and  samples  are  obtained  from  Equation  A.l  and  the  first  vk 
sample,  vq.  Even  with  this  initialization  there  is  a  transient  response  because  >o  and  jcq 
do  not  have  a  “history”  as  they  do  after  steady  state  is  achieved.  It  is  therefore 
suggested  that  the  filter  be  “warmed  up”  for  at  least  one  decorrelation  time  before 
using  the  output. 

An  filter  can  be  created  by  cascading  three  filters.  The  filter  equations 
are  therefore  given  by: 


zk  -  azk-i  +  by>j^_i 


yk  =  ayk-i  +  bxk-i  (A.l) 

xk  =  axk-i  +  bvk..]  . 


The  coefficienis  a  and  b  are: 

a  =  expf-tM^^^o)  =  expC-o^/No) 

where  =  2.904630.  Discrete  samples  of  the  additive  white  Gaussian  noise  random 
process  v*  are  generated  using  the  equation: 

L 


(l-gyPod-/?) 
l+4o^  +  a‘’  . 


Again  the  vk  samples  must  have  mean  power  given  by  the  quantity  in  the  square 
brackets  so  that  the  mean  power  of  the  filter  output  samples,  zh  will  have  mean  power 
of  The  discrete  channel  impulse  response  function  is 

h{kAt)  =  Zk  +  V^o^e*^  . 

Again,  to  minimize  the  transient  response  at  start-up  the  initial  filter  values  aie: 

Zo  =  yj  Pq{1—R)  uq 


>'0  = 


(l-aVo(l-/?) 


L 

2 


L  \+4a^  +  a*  J 


“l 


A~4 


Xo  = 


(l-aVPad-/?) 
\+Aa^  +  a* 


«2 


where  uq,  Ui,  and  U2  are  independent  samples  of  complex  AWGN  uncorrelated  with  the 
vk  samples.  The  filter  should  be  “warmed  up”  for  at  least  one  decorrelation  time 
before  the  output  is  used. 

A.3  WHAT  ABOLT/"^  PSDs? 

An  obvious  question  is:  Why  not  discuss  generation  of  realizations  with  an 
Doppler  PSD?  Indeed  the  expressions  for  xk  in  Equations  A.  I  and  A.2  are  single-pole 
filters  that  produce  realizations  with/"^  Doppler  PSDs: 

Xk  =  axk-i  +  bvk-i  . 

For  this  simple  case,  the  coefficients  a  and  b  are: 

a  =  expi-At/Zo)  =  exp(-l/iVo) 

b  =  . 

The  problem  with  this  Markov  process  is  that  the  temporal  statistics  are,  strictly 
speaking,  undefined.  This  is  because  the  scale  factor  A  in  the  expression  for  the  mean 
number  of  level  crossings  in  Section  2.4  is 

oo 

i  I  ito(ODfS{o)D)^  .  (A.3) 


For  the  Doppler  PSD, 


S(cod) 


_^Sl_ 

1  +  TlQ)j)  ' 


and  the  expression  for  A  yields  infinity. 

Fortunately,  for  a  discrete  realization  there  is  a  maximum  Doppler  frequency 
determined  by  the  sample  spacing.  The  highest  Doppler  frequency  component  in  the 
sampled  realizations  has  a  period  2At,  conesponding  to  a  maximum  frequency  of: 

(^D,max  ^ 


The  value  of  A  for  a  sampled  realization  is  then: 

tan'*  (tiA^o) 


A-  =  N„- 


n 


A~5 


For  integer  values  of  Nq  this  expression  reduces  to  Nq. 

Thus  the  mean  number  of  level  crossing  is  finite  and  the  mean  duration  and 
separation  of  fades  aie  non-zero.  Unfortunately  these  quantities  depend  on  Nq.  As  the 
number  of  samples  per  decorrelation  time  is  increased,  the  mean  number  of  level 
crossing  increases  and  the  mean  fade  duration  decreases. 

The  dependence  of  the  temporal  statistics  of  the/~^  Doppler  PSD  realization  on 
Nq  is  likely  to  be  unacceptable  in  most  applications.  Thus  only  Doppler  PSDs  with  a 
frequency  roll-off  greater  than  /"^  (so  Eqn.  A.3  is  finite)  should  be  used  to  generate 
realizations  of  Rician  fading.  This  will  ensure  that  the  realization  temporal  statistics  are 
well-behaved. 


APPENDIX  B 

JOINT  PROBABILITY  DENSITY  FUNCTION  /(a^  ) 


The  purpose  of  this  appendix  is  to  derive  the  joint  probability  density  function  of 
the  Rician  amplitude  a  and  its  time  derivative  a'  =  da/dt.  This  function  is  required  to 
calculate  the  temporal  statistics  of  Rician  fading.  A  less  general  form  of  this  derivation 
was  first  publish*^  by  Rice  [1948]. 

The  starting  point  for  this  calculation  is  the  determination  of  the  joint 
probability  density  function  of  the  random  in-phase  and  quadrature  components  x  and 
y  of  the  complex  envelope  of  the  electric  field.  It  is  assumed  that  x  and  y  are 
independent,  have  zero  mean,  and  that  they  are  normally  distributed.  Thus  the  joint 
probability  density  function  of  v  and  y  is 


fix,y)  = 


x^+\^ 
■  2(^ 


(B.l) 


Now  the  joint  probability  density  function  of  the  time  derivatives  x'  =  dx/dt  and  y'  ~ 
dy/dt  must  be  calculated.  It  will  be  shown  that  x  and  x'  are  independent,  as  arc  y  and 
y‘.  It  will  be  assumed  that  x,  x\  y,  and  y'  are  jointly  independent.  Thus  the  joint 
probability  density  function  of  x'  and  y  is  all  that  is  needed  in  addition  to  Equation  B.l 
to  write  down  the  joint  probability  density  function  J[xy,yy).  Once  this  function  has 
been  obtained,  a  simple  change  of  variables  from  x,  x',  y,  and  >•'  to  a  and  a'  will  yield 
the  desired  function. 


In  order  to  determine  the  distribution  of  x'  (or  y')*  consider  the  random  function 
jc(/)  written  as  a  Fourier  stochastic  integral 


xCr) 


oo 


•^oo 


(cod)  exp  (icoDt) 


dCDp 

2k  ■ 


(B.2) 


The  quantity  zio^n)  is  a  random  function  in  the  Doppler  frequency  domain.  It  is  useful 
to  assume  that  zicoo)  is  a  zero-mean,  normally  distributed  random  process,  although 
this  is  not  necessary  because  the  central  limit  theorem  will  make  xit)  normally 
distributed  for  almost  any  reasonable  disuibution  of  z{cop).  However,  with  the  normal 
assumption  for  z(Q)d),  x{t)  is  the  sum  of  many  independent,  normally  distributed 
random  variables,  and  is  necessarily  a  zero-mean,  normally  distributed  random 
variable. 

Before  continuing,  it  is  interesting  to  show  the  relationship  between  the  random 
spectral  components  z{Q)d)  and  the  Doppler  spectmm  Sicop).  The  correlation  function 
of  the  stationary  process  x(t)  may  be  written  as 


pih-h)  = 


<x(ri)jr(f2)) 


(B.3) 


dCiDn\ 

2n 


da)D2  {z(0}Di)Z*{Q)D7)) 

2n  ^2  QT^piicooih-icODzh)  • 


However  the  correlation  function  p(r)  may  also  be  written  in  terms  of  S{q)d): 


Pit)  = 


I  S(cod)  exp  iia)Dt)  ~  . 


(r‘.4) 


The  spectrum  Sicoo)  must  be  an  even  function  if  the  correlation  function  p(T)  is  to  be 
real. 


To  ensure  the  integral  in  Equation  B.3  is  only  a  function  of  time  difference 
t=ti-i2,  the  integrand  must  contain  a  factor  27c5( 6)01-0)02).  Using  the  Dirac  delta 
function  to  collapse  the  double  integral  in  Equation  B.3  and  comparing  the  result  with 
Equation  B.4  gives 

(zi(0Di)z*(C0D2))  =  2X0^  5(0)01-6)02)  5(0)oi)  .  (B.5) 

This  equation  also  demonsbates  that  the  random  Doppler  spectral  components  of  z(O)o) 
are  uncorrelated,  which  is  a  consequence  of  the  assumption  that  the  random  process 
x(r)  is  stationary. 

The  time  derivative  of  x{t)  is  given  by  differentiating  Equation  B.2,  with  a 
similar  expression  holding  for  y': 


At) 


00 

J (/6)o)  z(6)o)  exp(/0)oO 


dcop 
2k  ■ 


Because  z(q)o)  is  normally  distributed  with  zero  mean,  x'it)  will  also  be  normally 
distributed  with  zero  mean.  The  variance  of  x'it)  is 


{At)At)) 


00 


—  QC 


*wOO 


COD\(t^2iz.iQ}Di)z*iCOD2))  CXp  [i(0toi-0)o2)t] 


= 


7  dCJn 

cooSimo)^ 


=  -o^ 


d^pjx) 

d^ 


T=0 


(B.6) 
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The  variance  of  x'(0  may  be  written  in  the  general  fonn: 


where 


A  = 


1 

^  =  1.1858 

^tVV2  =  1.5176 


Gaussian  Doppler  PSD 
/~^  Doppler  PSD 
f~^  Doppler  PSD 


and  where  the  parameters  04  and  a«  were  determined  in  Section  2.2. 
The  cross  correlation  of  x'(t)  and  x(t)  is 


(x'(r)x(r)) 


V&)Dt 

2n 


(i6)oi) 


fdCDD2 

J  2k 


<—00 


(z(ahik*(o)D2))  U(o)Di~o)o2)r] 


= 


J  Ci)dS(COd) 


d(0[) 

2k 


=  0  . 


(B.7) 


Equation  B.5  and  the  fact  that  S(tt)o)  is  an  even  function  have  been  used  in  reducing 
Equations  B.6  and  B.7.  Because  x(r)  and  x'(0  are  uncorrelated  and  normally 
distributed,  they  are  also  independent.  Identical  results  hold  for  the  variance  of  y'  and 
the  cross  correlation  of  y  and  y'. 

The  joint  probability  density  function  of  x,  x',  y,  and  /  may  now  be  written 

down: 


{  ‘  1 

r 

(  -p-  y 

Poix^+y'^) 

exp 

20^  . 

4k(Pa^ 

exp 

4cPa^ 

f{xX,yy)  - 


This  function  may  be  U’ansform  to  the  desired  function  of  a  and  a  by  making  the 
change  of  variables 


x  +  r  cos  t?  =  a  cos  0 
y  +  r  sin  t?  =  a  sin  0 

where  r  is  the  constant  component  of  Rician  fading  and  t?  is  the  constant  phase.  The 
time  derivatives  of  x  and  y  are 


B-3 


x'  =  a'  cos  Q-a&  sin  6 
y'  =  a  sin  6  +  aff  cos  9 


which  gives  the  polar  coordinate  equations 

x^  +  y^  =  a^  +  -  2ar  cos  (6- 1?) 

The  probfaoility  density  function  coordinate  transformation  is 

jcy,y,y')  dxdx'dydy'  =  j{a,a',6^ff)  IdetfJ)!  dada'dGdff 
where  the  determinate  of  the  Jacobian  of  the  transformation  is 


det(i)  =  det 


"  ^ 

K 

d_i 

a/-' 

da 

da 

da 

da 

dx 

dx' 

da' 

do' 

da' 

da' 

dx 

djc' 

hi 

30 

30 

do 

30 

dx 

d/ 

-  dff 

y 

30' 

30'  J 

=  det 


= 


cos6  -ffsind  sin0  ffcosO 

0  COS0  0  sin0 

-asin0  -a'sin6-a0'cos0  acos9  acos9-affs\n6 
L  0  -asin0  0  ocos0 


The  joint  probability  density  function yfa,a', 0,0')  is 

~2ar  cos  (0- 1>) 


fya,a'Ad')  = 


Ind^ 


exp 


2a^ 


Tg(a'2  +  a0'2) ' 

4k^A^<^ 

\  > 

exp 

44V 

The  joint  probability  density  function  of  a  and  a'  is  obtained  by  integrating  this 
equation  over  0  and  ff : 


B-4 


with  the  result 


2jc  “> 

Jia,a')  =  jde  !d9'f{a,a',e,er)  , 

0  — oo 


yio.a') 


fa] 

r  + 

ar 

i?J 

exp 

20^  . 

f  T  y 

exp 

(B.8) 


lo  is  the  modified  Bessel  fuuction  that  results  from  performing  the  integral: 

Io(z)  =  Jexp(2Cos0)J0  . 

Thus  it  is  apparent  from  Equation  B.8  that  the  probability  density  function  of  a  is 
Rician;  the  probability  density  function  of  a'  is  normal  with  zero  mean  and  variance  of 
24^<7^/to;  and  a  and  a'  are  independent  because  theii'  joint  probability  density  function 
is  separable  into  a  function  of  a  times  a  function  of  a. 


B-5 


APPENDIX  C 

VARIATION  IN  THE  MEASUREMENT  OF  MEAN  POWER 

The  moment  of  the  amplitude  ak  of  an  impulse  response  function  realization 
with  N  samples  is  measured  using  the  formula 


Because  of  the  finite  length  of  a  realization,  pn  is  a  random  variable.  The  purpose  of 
this  appendix  is  to  develop  general  expressions  foi  the  mean  and  variance  of  P/j,  and 
then  to  apply  those  expressions  to  compute  the  expected  variation,  of,  in  the  measured 
mean  power  of  a  realization: 

Op  =  • 

This  variation  depends  primarily  on  the  number  of  decorrelation  times  in  the 
realization,  N/Nq,  where  Nq  is  the  number  of  samples  per  decorrelation  time.  It  is 
weakly  dependent  on  tbe  Doppler  frequency  power  spectral  density  (PSD)  and  on  the 
values  of  N  and  Nq. 

The  mean  value  of  pn  is  easy  to  compute: 

<Pn>  =  =  <a^) 

where  (a”)  is  the  ensemble  mean  value  of  the  moment  of  amplitude.  The  second 
moment  of  Pn  is  a  little  more  of  a  challenge  to  compute: 

1  V  <a^)  2 

<|i5>  =  <❖">  = 

where  Rrj(k)  is  the  correlation  of  the  moment  of  amplitude: 

For  the  general  case  of  Rician  fading  and  arbitrary  n,  the  joint  probability  density 
function  of  the  mplitude  at  two  times, /(a,, 02)*  is  needed  to  compute  the  correlation 
function.  However  in  the  special  case  of  mean  power  where  n  is  two,  the  correlation 
function  is  easily  computed  from  the  statisucs  of  the  underlying  complex  voltage. 

The  power  in  a  sample  of  the  impulse  response  function  with  a  Rician  amplitude 
distribution  is 


^  =  (^it  +  cos  +  [yjt  +  r  sin 


(C.l) 


C-I 


where  xk  and  yk  are  unconelated,  normally-distributed  random  processes  with  ^ero- 
mean  and  variance  The  “Rician”  components  r  cost?  and  r  sint>  are  constant.  It  is 
assumed  that  the  random  processes  x  at  two  times  is  jointly  normal: 


1 

- 

2jta"''l-p^ 


X^-lpX^Xi  +  xl 
2a\\-p'^) 


Values  of  the  two-point  corr^lfc  on  p  are  determined  by  the  functional  form  of  the 
Doppler  PSD  and  the  time  difference  between  the  samples  jtj  and  X2.  A  similar 
expression  holds  for  the  joint  probability  density  function  of  yi  and  y2. 

To  compute  the  variance  of  the  mean  power  measurement,  the  quantity  {a\(^  is 
required.  Using  Equation  C.l  this  quantity  involves  terms  of  the  form 

<Xi>  =  <X2>  =  <>!>  =  <>2)  =  0 

=  (-^i)  =  <vi>  =  (yl)  =  cr^ 


ixxxi)  =  (yiya)  = 

=  {x\x^  =  {y,y^)  =  (yfyj)  =  0 

<44)  =  (yhl)  =  +  2p^)  . 

The  CTOSS  correlation  <aia|)  is  then 

=  4o'‘(l  +p^)  +  4r^o^(l  +  p)  +  /  . 

The  Rician  amplitude  r  and  the  variance  are  written  in  terms  of  the 
scintillation  index  S4  so  that  the  mean  power  (a'‘)  is  constant  and  equal  to  Pq. 

=  PqR 

2a^  =  Poil-R) 


where  the  “Rician”  index  is 


R  = 


Combining  these  results  gives  the  following  expression  for  the  variance  in  the 
measured  mean  power: 


,  N-\ 

i-NR^  2  y 
N 


{\-k/N)[R  +  {\-R)p{k)f  . 


(C.2) 


The  two-point  correlation  p{k)  is 


C-2 


exp  [-(fcWo)^  ]  Gaussian  Doppler  PSD 

p(jfc)  =  [l  +  +  ]exp  [-OfiWo]  /“^  Doppler  PSD 

,  [1  +  a^kyNQ\Q\^  {-ank/No]  f~^  Doppler  PSD 

where  the  coefficients  04  (a4=2.146139)  and  (06=2.904630)  are  determined  by  the 
requirement  that  p(//o)=6"'- 

Plots  of  the  power  measurement  standard  deviation  are  shown  in  Figure  20  for 
N  equal  to  1024,  2048,  and  4096.  The  value  of  Nq  is  10  for  each  case,  so  the  curves 
conespond  to  realizations  of  length  102.4%  204.8%  and  409.6%  respectively.  Solid 
lines  in  the  figure  are  for  a  Gaussian  Doppler  frequency  PSD,  and  solid  circles  are  for 

Doppler  PSD. 

As  expected,  the  mean  power  measurement  variation  is  larger  for  realizations 
with  fewer  decorrelation  times.  The  measurement  variation  decreases  with  decreasing 
scintillation  index  because  the  fluctuating  part  of  the  impulse  response  contributes  less 
aud  less  to  the  total  power.  For  a  given  value  of  N  there  is  little  difference  between  the 
results  for  the  two  PSDs.  This  is  because  the  three  correlation  functions  above  differ 
little  for  values  of  k  less  than  Nq  where  p(A:)  is  close  to  unity,  but  vary  significantly  for 
larger  values  of  k  where  p()i:)  is  small.  Thus  the  differences  in  p()t)  for  differing 
Doppler  frequency  PSDs  occur  in  Equation  C.2  at  values  of  k  that  couuibute  little  to 
the  sum. 
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j  Scintillation  Index,  S> 
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j  Figure  20,  Power  measurement  error. 
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ATTN:  NTB/JPO  DR  C  GIESE 

PHILLIPS  LABORATORY 
ATTN:  NTN 

UNITED  STATES  STRATEGIC  COMMAND 
ATTN:  J  533 
ATTN:  J  534 
ATTN:  J  614 

DEPARTMENT  OF  ENERGY 

EG&G,  INC 

ATTN:  D  WRIGHT 

UWRENCE  LIVERMORE  NATIONAL  LAB 
ATTN:  L-97TOONICH 

LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  DANWINSKE 
ATTN:  RW  WHITAKER 

SANDIA  NATIONAL  LABORATORIES 
ATTN:  D  DAHLGREN  6410 
ATTN:  DIV  2344  ROBERT  M  AXLINE 
ATTN:  ORG  1231  J  R  LEE 
ATTN  ORG  9110  G  CABLE 
ATTN:  ORG  9110  WD  BROWN 
ATTN:  TECH  LIB  3141 

OTHER  GOVERNMENl 

central  INTELLIGENCE  AGENCY 
ATTN:  OSWR/NED 
A1TN:  OSWR/SSD  L  BERG 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORP 

ATTN:  BRIAN  PURCELL 
ATTN:  C  CREWS 
ATTN:  C  RICE 
AITN:  DR  J  M  STRAUS 
ATTN:  G  LIGHT 
ATTN;  J  THACKER 
ATTN;  MROLEN2 

AUSTIN  RESEARCH  ASSOCIATES 
ATTN:  n THOMPSON 


AUTOMETRIC,  INC 

ATTN:  C  LUCAS 

BDM  INTERNATIONAL  INC 

ATTN.  W  LARRY  JOHNSON 

BERKELEY  RSCH  ASSOCIATES,  INC 
ATTN:  J  WORKMAN 
ATTN:  NTGLADD 
ATTN:  S  BRECHT 

DELVIN  SYSTEMS 

ATTN:  BPHll.LIPS 
ATTN:  NCIANOS 

DYNETICS,  INC 

ATTN:  WILLIAM  D  TEPPER 

ELECTROSPACE  SYSTEMS,  INC 
ATTN:  LINDA  CALDWELL 
ATTN:  P  PHILLIPS 

EOS  TECHNOLOGIES,  INC 
ATTN:  B  GABBARD 
ATTN;  RLELEVIER 

GENERAL  RESEARCH  CORP  INC 
ATTN;  J  EOLL 

GRUMMAN  AEROSPACE  CORP 
ATTN:  J  DIGLIO 

H.ARRIS  CORPORATION 
ATTN:  EKNICK 
ATTN;  LYMUEL  MCRAE 

HORIZONS  TECHNOLOGY,  INC 
ATTN;  BLEE 

INFORMATION  SCIENCE,  INC 
ATTN:  WDUDZIAK 

INSTITUTE  FOR  DEFENSE  ANALYSES 
ATTN-  E  BAUER 
ATTN:  HWOLFHARD 

JAYCOR 

ATTN:  J  SPERLING 

KAMAN  SCIENCES  CORF 
ATTN;  DASIAC 
ATIN;  E  CONRAD 
ATTN:  GDITTBERNER 

KAMAN  SCIENCES  CORPORATION 
ATTN:  BGAMBILL 
ATTN:  DASIAC 
ATTN:  R  RUTHERFORD 

LOCKHEED  MISSILES  &  SPACE  CO,  INC 
ATTN:  J  KUMER 
ATTN:  R  SEARS 

LOGICON  R  &  D  ASSOCIATES 
ATTN.  D  CARLSON 
ATTN:  S  WOODFORD 


Di3t-2 


DNATR-92'98  (DL  CONTINUED) 


LOGICON  R  &  D  ASSOCIATES 
ATTN;  J  WALTOf; 

LOGICON  R  &  D  ASSOCIATES 
ATTN:  EHOYT 

MAXWELL  LADS.  INC 

ATTN:  BILLRIX 

MCDONNELL  DOUGLAS  CORP 
ATTN:  TCRANOR 

MCDONNELL  DOUGLAS  CORPORATION 
ATTN:  RHALPRIN 

MISSION  RESEARCH  CORP 
ATTN;  J  KENNEALY 
ATTN:  R  ARMSTRONG 
ATTN:  S  DOWNER 
ATTN:  W  WHITE 

MISSION  RESEARCH  CORP 
ATTN:  RLBOGUSCH 

MISSION  RESEARCH  CORP 
ATTfJ;  DAVEGUICE 

MISSION  RESEARCH  CORP 
ATTN:  BR  MILNER 
ATTN:  B  SAWYER 
ATTN;  DKNEPP 
ATTN:  DUNDMAN 
A^TTTiI-  C  GUIGLIANO 
ATTN:  RBIGONI 
2  CYS  ATTN:  R  DANA 

ATTN;  R  HENDRICK 
ATTN:  SGUTSCHE 
ATTN:  TECH  LIBRARY 

MITRE  CORPORATION 

ATTN:  DR  MR  DRESP 

MITRE  CORPORATION 

ATTN:  G  COMPARETTO 

NORTHWEST  RESEARCH  ASSOC.  INC 
ATTN:  E  FREMOUW 

PACIFIC-SIERRA  RESEARCH  CORP 
ATTN:  H  BRODE 
ATTN:  RLUTOMIRSK.I 

PACIFIC-SIERRA  RESEARCH  CORP 
ATTTJ:  MALLERDING 


PHOTOMETRICS,  INC 

AITN:  I  L  KOFSKY 

PHOTON  RESEARCH  ASSOCIATES 
ATng;  DBURWELL 

RJO  ENTERPRISES/POET  FAC 

ATTN;  STEVEN  KRAMER 

S-CUBED 

ATTN;  C  NEEDHAM 

SCIENCE  APPLICATIONS  INTL  CORP 
«TTN:  C  SMITH 
ATTN:  D SACHS 
2  CYS  ATTN:  L  LINSON 

SCIENCE  APPLICATIONS  INTL  CORP 
A ITN:  J  COCKAYNE 

SPARTA  INC 

ATTN:  KCOSNER 

SPARTA  INC 

ATTN.  DDEAN 

SR!  INTERNATIONAL 

ATTN:  R  LIVINGSTON 
ATTN:  RTTSUNODA 
ATTN:  WCHESNUT 

STEWART  RADIANCE  UBORATORY 

A-TTKJ.  p  H.UPPI 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  J  FORD 
ATTN:  J  VVOLFSBERGER  JR 
ATFN:  N.'ASSINO 
ATTN:  RONALD  E  LEWIS 

TOYON  RESEARCH  CORP 
ATTN;  J  ISE 

TRW  SPACE  &  DEFENSE  SECTOR 
ATTN;  DM  LAYTON 
ATTN:  HL  DEPT  LIBRARY 

VISIDYNE,  INC 

ATTN:  J  CARPENTER 
ATTN-  J  DEVORE 
ATTN:  J  T.HOMPSON 
ATTN;  W5CHLUETER 
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